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We give a pedagogical introduction to two aspects of magnetic fields in the early universe. We first focus on how to 
formulate electrodynamics in curved space time, defining appropriate magnetic and electric fields and writing Maxwell 
equations in terms of these fields. We then specialize to the case of magnetohydrodynamics in the expanding universe. We 
emphasize the usefulness of tetrads in this context. We then review the generation of magnetic fields during the inflationary 
era, deriving in detail the predicted magnetic and electric spectra for some models. We discuss potential problems arising 
from back reaction effects and from the large variation of the coupling constants required for such field generation. 
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1 Introduction 

The universe is magnetized, from scales of planets to the 
largest collapsed objects like galaxy clusters. The origin and 
evolution of these magnetic fields is a subject of intense 
study. Much of the work on magnetic field origin centers 
on the idea that small seed magnetic fields are generated by 
purely astrophysical batteries and are subsequently ampli- 
fied to the observed levels by a dynamo (cf. Brandenburg 
and Subramanian, 2005 for a review). An interesting alter- 
native is that the observed large-scale magnetic fields are 
partially a relic field from the early Universe, which have 
been further amplified by motions. 

We provide here a pedagogical review of two aspects of 
magnetohydrodynamics (MHD) in the early universe. In the 
first part of the review we focus on how to formulate elec- 
trodynamics in curved space time, especially how to define 
magnetic and electric fields and write Maxwell equations in 
terms of these fields. This issue may perhaps be well known 
to relativists but may not be so familiar to practitioners of 
MHD. The perspective provided is that of the author. 

In the second part, we review, again in a pedagogical 
manner, generation of magnetic fields during the inflation- 
ary era. It is well known that scalar (density or potential) 
perturbations and gravitational waves (or tensor perturba- 
tions) can be generated during inflation. Could magnetic 
field perturbations also be generated? Indeed, a large num- 
ber of mechanisms whereby magnetic fields are generated 
in the early universe have been discussed in the literature 
(Turner & Widrow 1988; Ratra 1992; Widrow 2002; Gio- 
vannini 2007). These generically involve the breaking of 
the conformal invariance of the electromagnetic action, and 
the predicted amplitudes are rather model dependent. Nev- 
ertheless, if a primordial magnetic field is generated, with a 



present-day strength of B ~ 10 -9 G and coherent on Mpc 
scales, it can strongly influence early galaxy formation (cf. 
Sethi and Subramanian, 2005), or induce signals on the cos- 
mic microwave backround radiation (CMBR) (cf. Subrama- 
nian, 2006 for a review) including CMBR non-Gaussianity 
(cf. Seshadri and Subramanian, 2009; Caprini et al., 2009). 
An even weaker field, sheared and amplified due to flux 
freezing, during galaxy and cluster formation may kick-start 
the dynamo. It is then worth considering if one can generate 
such primordial fields. Our discussion of inflationary gener- 
ation of magnetic fields follows the standard literature, but 
again the perspective offered is that of the author. We hope 
that our pedagogical discussion of these aspects of MHD in 
the early universe will be useful to those entering the sub- 
ject. 

2 Electrodynamics in curved spacetime 

We discuss to begin with Maxwell equations in a general 
curved spacetime and then focus on FRW models. Elec- 
trodynamics in curved spacetime is most conveniently for- 
mulated by giving the action for electromagnetic fields and 
their interaction with charged particles: 

1 



S = - V=g 



Here F = A 



— g d x 



(1) 



= A VtfJl — A^ v is the electromag- 
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netic (EM) field tensor, with A^ being the standard electro- 
magnetic 4-potential and J M the 4-current density. Further, 
here and below, we use Greek indices /i, v etc for spacetime 
co-ordinates and Roman indices i,j, k... for purely spatial 
co-ordinates; repeated indices are summed over all the co- 
ordinates. We also adopt units where the speed of light c = 
1 and a metric signature (—,+,+,+). Demanding that the 
action is stationary under the variation of A^, gives one half 
of the Maxwell equations, 

F» v . v = 4irJ^ (2) 
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And from the definition of the electromagnetic field tensor 
we also get the other half of the Maxwell equations 

F^u; 7] = 7] = 0; ov*F^. l/ = 0. (3) 
The square brackets \\iv 1 7] means adding terms with cyclic 
permutations of /i, v, 7. In the latter half of Eq. (O, we have 
defined the dual electromagnetic field tensor 



1 iivuPp 

2 ap 



Here ^ va is the totally antisymmetric Levi-Civita tensor, 
1 \ ' _ — , 

£pvp\ — V Q-^-pvpXi 



vp\ 



and g is the determinant of the metric tensor. Further A^ vpX 
is the totally antisymmetric symbol such that .4 0123 = 1 
and ±1 for any even or odd permutations of (0,1,2,3) re- 
spectively. Note that -4oi23 = — 1- 

We would like to cast these equations in terms of elec- 
tric and magnetic fields (Ellis, 1973, Tsagas, 2005, Barrow, 
Maartens & Tsagas, 2007). In flat spacetime the electric and 
magnetic fields are written in terms of different components 
of the EM tensor F„„. This tensor is antisymmetric, thus 
its diagonal components are zero and it has 6 independent 
components, which can be thought of the 3 components of 
the electric field and the 3 components of the magnetic field. 
The electric field E l is given by time-space components of 
the EM tensor, while the magnetic field B l is given by the 
space-space components 

F 0i = pi pl2 = B 3 p23 = B l p31 = B 2 (4) 

In a general spacetime, to define corresponding electric and 
magnetic fields from the EM tensor, one needs to isolate a 
time direction. This can be done by using a family of ob- 
servers who measure the EM fields and whose four-velocity 
is described by the 4-vector 

vr = — — ; uru,, = — 1 . 
ds ^ 

Given this 4-velocity field, one can also define the 'projec- 
tion tensor' 

This projects all quantities into the 3-space orthogonal to u M 
and is also the effective spatial metric for these observers, 
i.e 

ds 2 = g^dx^dx" = -(u^dx^ 1 ) 2 + h^dx^dx" . (5) 

Using the four-velocity of these observers, the EM fields 
can be expressed in a more compact form as a four-vector 
electric field and magnetic field as 

1 



Ep = F^i 



From the definition of E^ and B^ 



Epu» 



»F" x = *F^u v . (6) 
we have 

0, B^u" = , 
Thus the four-vectors £?„ and E^ have purely spatial com- 
ponents and are effectively 3-vectors in the space orthogo- 
nal to u p . One can also invert Eq. (0 to write the EM tensor 
and its dual in terms of the electric and magnetic fields 



F, 



UpE v 



u„E„ 



B a vP 



(7) 



* pafi 



_ e aPp,i> p 

2 



— c appv 



u„E v + {^B? - B a vP) . 



(8) 



We can now use the time-like vector u M and the spatial met- 
ric h% to decompose the Maxwell equations into timelike 
and spacelike parts. Consider projection of Eq. (0 on u a . 
We have 



u a (*F a0 .p) =0 = (u a *F a ^),p - u a , p * F a ^ 



(9) 



Substituting Eq. ® into Eq. (O we get 

[u^^u^Ev + u a u a B$ - u u a B a ] 

-u a -p [e a ^ v UfJ ,E v + u a BP - u^B a ] = (10) 

The 1st, 3rd and 5th terms are zero because respectively 
Uae^^Uf, = 0, u a B a = and u a . p u a = 0. The rest 
of the terms give 



B 







(11) 



where we have defined the acceleration 4-vector, 



Up =u up- a 

as the directional derivative of u 13 , along the timelike direc- 
tion specified by u a . The covariant velocity gradient ten- 
sor, u a .p, can be decomposed into, shear, expansion, vortic- 
ity and acceleration parts in the following manner: We first 
write 

u a -p = 8£up- tl/ [5p + u v up\ - 5^u l _ l . u u u u f} 
= h£hpU KI/ - u a up 

= ©a/? + W Q/3 - U a Ufj (12) 

In the last line, we have decomposed the spatial part of the 
co-variant derivative h^hpU^, into symmetric expansion 
tensor <d a p and an antisymmetric vorticity tensor uj a p. Note 
that both the expansion and vorticity tensors are purely spa- 
tial, in the sense that <d a puP — — u> a pu^ . One can further 
split Q a p into its trace and trace free part. 

Qa/3 = CT a p + ~@h a p 

where 9 = 0„ = u a a is called the expansion scalar, and 
o a p is called the shear tensor. Note that the shear tensor 
satisfies, a" = 0, and a a pu° = 0, that is it is traceless and 
also purely spatial. Thus we have 



u a -p = a a p + -@h a p + u a p — ii a up 



(13) 



Only the antisymmetric part of u a -p contributes when we 
substitute Eq. ( fT3l into Eq. ( fTTT >- Further simplification can 
be made by defining the vorticity vector, 

u v = -\u a pe a ^ Ufi 
and the spatial projection of the covariant derivative 

DpB a = h%KB^ . 
Then we have 

DpB* = h%hlB^ = KBl = (8S + u»u v )Bl 

= B£ - vTu^B" . (14) 
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Thus Eq. ( fTTT i becomes 

DpB fj = 2w /3 Ep . (15) 

This equation generalizes the flat space equation V • B = 0, 
to a general curved spacetime. We see that 2ui^Ep acts as 
an effective magnetic charge, driven by the vorticity of the 
relative motion of the observers measuring the electromag- 
netic field. 

Now turn to the spatial projection of of Eq. ([3]l on h*. 
We have 

h%{*F aP . p ) = h K a [e a ^ Ull E v +u a Bf> - U PB a ]. 
= (16) 
where we have again substituted Eq. ^ for the dual EM 
tensor. Expanding out the covariant derivative in Eq. dT6b 
we get 



"a 

hi 



z a ^(u K pE v 



B a 



u^E v ,p) + u^B* 







u a B. 



(17) 



Using h* = (6% + u K u a ) above, only the contribution of 



8* survives in the 2nd term (since e 



ap\iv 



0), 3rd 



term (since u a u a ^ = 0) and the 5th term (since u a B a = 0), 
while the 4th vanishes because h*u a = 0. Thus we have 
from the remaining terms, 

Ke a ^ v u^E v + e K ^" UlM E v;0 



-,0 



B? - u'^B K - h*B a = 0(18) 

The first term can be simplified using Eq. ( fT3l . Due to the 
antisymmetry of the Levi-Cevita tensor, the symmetric parts 
in Eq. dT3T > drop out and we have, 

K^u^E v = h K a e a ^"E^p 

- h K a € a ^ v u^upE u . (19) 

The first term on the RHS of the above equation vanishes, i.e 
^ v E v u) up — 0. To see this it is convenient to define a 
3-d fully antisymmetric tensor 



c kPv/j, 



Note that this tensor is also purely spatial in the sense that 
e K P u u v = I K ^ u up = e Kpy u K = 0, a result which follows 
from the anti-symmetry of We can then write the 4-d 

Levi-Civita tensor as 

(20) 

A 0a )/2. The 



where we use the notation, A^ a ^ ■■ 
first term in Eq. ( fT9l l then becomes 

hge^E^p = h K a E v u^ x 



(-4 



a/3 



we 



" + l a ^u u ] = 



The last equality follows from the fact that h^u a 
E u v — n anA '■' -"'P 



(21) 
0, 



and uj^pu 1 * 



n„e 1 h u,,,rE, 



0.Q Therefore we have 

K^u^Ev 
e K ^ v u^upE v . (22) 



1 A more physical way of seeing that the first term in Eq. 1 1 9i vanishes 
is to note that E v and u)^p are purely spatial tensors, orthogonal to the 
'time' direction specified by u a . Thus only the time like component of 
pa _ e a Pi J - 1 ' EvUi^p is non-zero. This implies however that h^P a = 0, 
thus proving the required result. 



The second term of Eq. ( TT8l can be written more transpar- 
ently by defining a 'Curl' operator 

Curl(£ K ) = l Kliv E v .fi 

Thus e^^u^E^p = -Curl(E K ). The third and fourth 
terms in Eq. ( fTFt can also combined and rewritten using 
Eq. ( fT3l to give, 

% - Q6$] B? 
.2. 



u%B« 



B 



(23) 



Note that the acceleration term in Eq. ( fT3l does not con- 
tribute above because upB^ = 0. Putting all the above re- 
sults together, Eq. ( TT8l gives the generalization of Faraday 
law to curved spacetime, 

. 2. 



hlB a = 



B 13 



Curl(E K ) . 



(24) 



The other two Maxwell equations, involving source terms, 
can be derived in very similar manner. Note that if we map 
E — > —B, and B — > E, then the dual EM tensor is mapped 
to the EM tensor, that is *F fiV F^ u . We can use this 
symmetry between the EM tensor and its dual to read off 
the Maxwell equations involving the source terms from the 
source free equations, Eq. ( fT5l l and Eq. (124-b . We also note 
that in deriving Eq. ( [TBI ) and Eq. ( l24b . we changed the sign 
of all the terms appearing in Eq. (fTOt and Eq. (fT7t . Thus 
mapping E -» —B, and B — > E in Eq. ( fTBT l and Eq. d24] > 
respectively, and also changing the sign of the source term 
47rJ M -> —AnJ^, the Maxwell equations F^ v . v = AttJ^, 
in terms of the E^ and B^ fields, become 

DpE& = inp q - 2^ Bp , (25) 



hlE a 



l^p 



Curl(B K ) - Airf 



(26) 



Here we have defined the charge and 3-current densities as 
perceived by the observer with 4-velocity u a by projecting 
the 4-current density J M , along u a and orthogonal to u a . 
That is 



Pq = ~J^U p , J" 

Note that j K u K = 0. To do MHD in the expanding universe, 
we also need the relativistic generalization to Ohm's law. 
This is given by 

aF afl wp, or 
J a = P(f) ,w a + aE? f) . (27) 



h a TP — 



Here the symbol (/) stands for a fluid variable, that is w c 
is the mean 4-velocity of the fluid, h\ 



(/)/3 



(5| + w a wp) 

and E c ^^ = F a ^wp is the electric field as measured in the 
fluid rest frame. Also P(f) q and a are the fluid charge den- 
sities and conductivity as measured in its rest frame. Note 
that the fluid 4-velocity w a , need not be the 4-velocity u a of 
the family of fundamental observers used to define the EM 
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fields in Maxwell equations; indeed the conducting fluid 
will in general have a peculiar velocity in the rest frame of 
the fundamental observers. 

Further discussion on electrodynamics in curved space- 
time (using the 3+1 formalism) and how the different parts 
of the spacetime geometry affect the EM field can be found 
in Tsagas (2005). We now specialize to case of the expand- 
ing universe. 

2.1 Electrodynamics in the expanding universe 

Let us now consider Maxwell equations for the particular 
case of the expanding universe, with the metric that of spa- 
tially flat Friedmann-Robertson- Walker (FRW) spacetime, 

ds 2 = -dt 2 + a 2 {t) [dx 2 + dy 2 + dz 2 ] . (28) 

Here t is the proper time as measured by the fundamental 
observers of the FRW universe, while x, y, z are co-moving 
spatial co-ordinates. The expansion of the universe is de- 
termined by the scale factor a(t), and H{t) = a/a is the 
Hubble expansion rate (we have also defined d = da/dt). 
We choose u a corresponding to the fundamental observers 
of the FRW spacetime, that is u a = (1, 0, 0, 0). For such a 
choice and in the FRW spacetime, we have 

ii a = 0, w Q/3 =0, cr Q/3 =0, 6 = 3-. (29) 

Further, we can simplify h*B a as follows: 

KB a = (6% + u K u a )u^B« 

= u^B* + m k u 7 [(m q B q ) ;7 - u an ] 
= U^B* (30) 
Thus the Maxwell equations reduce to, 
B% = 0, E% = 47rp g , 

v?B^ + \QB K = -Curl(£ K ), 

v?E^ + \QE* = Curl(B K ) - Airf . (31) 

In the spatially flat FRW metric the connection co-efficients 
take the form 

Too = = r 0j =Tj k , T° = Sijaa, 

<C- (32) 

Using these Eq. (f3TT > can be further simplified as follows: 

dB ' - dE i _ ^ 
dx 1 dx % q ' 

-11 = I e * m M - 47r f . (33) 

a 3 dt 1 1 a Um dx m 

Here we have defined the 3-d fully antisymmetric symbol 
£ ijk> w i tn e i23 = !• These equations resemble the flat space- 
time Maxwell equations except for the presence of the scale 
factor a(t). 



The electric and magnetic field 4-vectors we have used 
above are referred to a co-ordinate basis, where the space- 
time metric if of the FRW form. They have the following 
curious property. Consider for example the case when the 
plasma in the universe has no peculiar velocity, that is w a — 
u a , and also highly conducting with a — s- oo. Then from 
Eq. ( |27] |. we have E"^ = = E a , and from Faraday's 
law in Eq. d33t . B l oc 1/a 3 . There is however a simple re- 
sult derivable in flat space time that in a highly conducting 
fluid, the magnetic flux through a surface which co-moves 
with the fluid is constant. Since in the expanding universe all 
proper surface areas increase as a 2 (t), one would expect the 
strength of a 'proper' magnetic field to go down with expan- 
sion as 1/a 2 . This naively seems to be at variance with the 
fact that B % oc 1/a 3 and Bi = gi^B^ oc 1/a. There are two 
comments to made at this stage: First, if we define the mag- 
netic field amplitude, say B, by looking at the norm of the 
four vector B^, that is let B 2 = B^B^ = B t B l oc 1/a 4 , 
then we do get B oc 1/a 2 . This procedure however does not 
appear completely satisfactory as one would prefer to deal 
with the field components themselves. Another possibility 
is to refer all tensor quantities to a set of orthonormal basis 
vectors, referred to as tetrads. 

Any observer can be thought to be carrying along her/his 
world line a set of four orthonormal vectors e( ), where a = 
0, 1, 2, 3, which satisfy the relation 

<Wef„)4)=^> r l abe Xa) e h=9^ (34) 
Here 77^ has the form of the flat space-time metric. We 
choose the observer's 4-velocity itself to be the tetrad with 
a — 0, i.e e^ Q % = it* 4 . The other three tetrads are orthogonal 
to the observer's 4-velocity. In the present case, we consider 
the observer to be the fundamental observer of the FRW 
space time, and the components of the tetrads, which satisfy 
Eq. d34T > are given by 

e (0) — °0> (i) a * ' 

The metric 77^ can also be used to raise and lower the index 
of the tetrad to define = rf'efty Note that the fun- 

damental observers move along geodesies, and as we noted 
earlier, do not have either relative acceleration or rotation. 
Such observers parallel transport their tetrad along their tra- 
jectory, i.e u^e? y = 0, as can be easily checked dy di- 
rect calculation using the connection co-efficients given in 
Eq. (l32t . The magnetic and electric field components can 
now be represented as its projection along the four orthonor- 
mal tetrads using, 

B a - g^B^W , E a = g^EWM , (35) 
which gives 

B° = 0, E° = 0, 

B a = a(t)B a , E a = a(t)B a , for(a = 1, 2, 3). (36) 

Note that B a , E a are co-ordinate scalars, but the set of four 
scalars B a is still a vector as far as local Lorentz transforma- 
tion is concerned (which preserves the orthonormality con- 
ditions in Eq. d34"i>). If we define B a = r) a i,B b , then numeri- 
cally B l = Bi and B° — — B = 0. Similar relations obtain 
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for the electric field components. In the FRW universe, as 
B % oc 1/a 3 , we see that B l = B t ex 1/a 2 , as one naively 
expects from flux freezing of the magnetic field. Thus the 
magnetic field components projected onto the orthonormal 
tetrads seem to be the natural quantities to be used as the 
'physical' components of the magnetic field. Note that this 
is similar to using the Cartesian components of a vector as 
the physical components in 3 dimensional vector analysis. 

There is an additional feature of using tetrads which is 
of particular interest. Given the set of tetrads one can set 
up a local co-ordinate system around any event V by us- 
ing geodesies emanating from V and whose tangent vectors 
at V are the four tetrads e( ). This co-ordinate frame, is a 
locally inertial frame; that is the spacetime is locally flat 
with the metric in the form of r) a b and the connection co- 
efficients in these co-ordinates vanishing (see section 13.6 
of Misner, Thome & Wheeler, 1971 (MTW) for a proof). 
In fact such a co-ordinate system can be set up all along 
the world line of the fundamental observer, and are called 
then Fermi-Normal co-ordinates (Manasse & Misner, 1963; 
MTW; Cooperstock, Faraoni and Vollich, 1998). To lead- 
ing order, the time direction in this inertial frame, is the 
proper time co-ordinate t of the FRW metric and the space 
co-ordinates become the proper space co-ordinates r % — 
a(t)x l . (There are second order deviations from these re- 
lations (cf. Cooperstock, Faraoni and Vollich, 1998)). 

The Maxwell equations in such a locally inertial frame 
take a particularly transparent form. They can be derived 
from Eq. OTb by adopting the flat space metric and replac- 
ing the covariant derivative with an ordinary derivative. As 
the metric is locally flat, the curl operator also reduces to 
the ordinary curl with respect to r % . Importantly 9 being a 
scalar is invariant under co-ordinate transformations, and so 
still = 3a/ a. The field components are the same as those 
defined above using the tetrads. This is because the tetrads 
become the co-ordinate basis vectors in the Fermi-Normal 
co-ordinates. The Maxwell equations then become, 



dB l 
Or 
0_ 
dt 
8 



= 0, 



d{a 2 E l ) 



Aitp q a , 



[a 2 B l ] = -e 



d{a 2 E l ) 



a 2 E i 



ilm Q r rn 
2 ul\ 



(37) 



In the absence of charges and currents, Eq. (1371 has elec- 
tromagnetic wave solutions, with the amplitude of the elec- 
tric and magnetic fields decaying with expansion as 1/a 2 . 
In the presence of a conducting medium, one has to again 
supplement these Maxwell equations with the Ohm's law. 
In the limit of non-relativistic fluid velocity v, this again 
reduces to 

f = PqV * + a[E* + e* Um Vl B m ] (38) 
If we neglect the displacement current and charge density 
terms, as valid for a highly conducting medium (cf. Bran- 
denburg and Subramanian, 2005), the induction equation 

becomes 
d 



where we have defined the vector B = (B 1 , B 2 ,B 3 ). Thus 
we see that in the absence of resistivity (77 = 0) or peculiar 
velocities (v = 0), the magnetic field defined in the local 
inertial frame, decays with expansion factor as B oc 1/a 2 . 
As pointed out above, this decays is as expected, when the 
magnetic flux is frozen to the plasma, since all proper areas 
in the FRW spacetime increase with expansion as a 2 . This 
completes our pedagogical discussion of doing magnetohy- 
drodynamics in curved spacetime and in particular the ex- 
panding universe. We turn to a pedagogical discussion of 
primordial magnetic field generation during the inflationary 
era. 

3 Magnetic field generation during inflation 

The early universe is supposed to have gone through an 
epoch of accelerated expansion referred to as inflation. In- 
flation can provide a solution to several problems of stan- 
dard big bang cosmology, one of them being to explain the 
origin of perturbations which eventually led to all the struc- 
tures that we see. We refer the reader to several standard 
textbooks for a discussion of the inflationary paradigm (Lid- 
die and Lyth, 1999; Mukhanov, 2005; Padmanabhan, 2002) 
and recent reviews (Bassett, Tsujikawa, Wands, 2006; Sri- 
ramkumar, 2009). Inflation provides several ideal conditions 
for the generation of primordial fields with large coherence 
scales (Turner & Widrow, 1988). First the rapid expansion 
in the inflationary era provides the kinematical means to 
produce fields correlated on very large scales by just the 
exponential stretching of wave modes. Also vacuum fluctu- 
ations of the electromagnetic (or more correctly the hyper- 
magnetic) field can be excited while a mode is within the 
Hubble radius and could be transformed to classical fluctua- 
tions as it transits outside the Hubble radius. Finally, during 
inflation any existing charged particle densities are diluted 
drastically by the expansion, so that the universe is not a 
good conductor; thus magnetic flux conservation then does 
not exclude field generation from a zero field. Most of the 
models for magnetic field generation during inflation take 
the field to be described by the action of an abelian gauge 
field, and have not considered the action obtained from for 
example the Electro- Weak or some Grand Unified theory. 
For simplicity we shall also adopt this approach below. 

There is one major difficulty, which arises when one 
considers magnetic field generation during inflation. This is 
because the standard electromagnetic action is conformally 
invariant, and the universe metric is conformally flat. Con- 
sider again the electromagnetic action 



S = - 



1 

167T" 

1 



g d*x —g^g^F, u F al3 



(40) 



Suppose we make a conformal transformation of the metric 
given by 



dt 



Be 



= V r x[t)X {Ba 2 ) - 77V x (Ba 2 )] (39) 



n 2 



9^v 



(41) 
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f=g* = 4 \/=5 and g*^ a = fT V"- Then 



This implies 
taking 

A*^ = A^S* = S. (42) 

Thus the action of the free electromagnetic field is invariant 
under conformal transformations. Note that the FRW mod- 
els are conformally flat; that is the FRW metric can be writ- 
ten as g»„ — £l 2 7]^ v , where rj^v is the Minkowski, flat 
space metric. As we show below explicitly for a universe 
with flat spatial sections, this implies that one can trans- 
form the electromagnetic wave equation into its flat space 
version. It turns out that one cannot then amplify electro- 
magnetic wave fluctuations in such a FRW universe and the 
field then always decreases @ with expansion as l/a 2 (t). 

Therefore mechanisms for magnetic field generation re- 
quire the breaking of conformal invariance of the electro- 
magnetic action, which changes the above behaviour to B ~ 
1 /a c with typically e< 1 for getting a strong field. A mul- 
titude of ways have been considered for breaking conformal 
invariance of the EM action during inflation. Some of them 
are illustrated in the action below: 

S = J^g^x [-ffafyJL-F^Ff-bRA 3 

+ gOF^F^ - D^{D»i>y ] (43) 

They include coupling of EM action to scalar fields ((f)) like 
the inflaton and the dilaton, coupling to curvature invariants 
(R), coupling to a pseudo-scalar field like the axion (9), hav- 
ing charged scalar fields (iji) and so on. If conformal invari- 
ance of the EM action can indeed be broken, the EM wave 
can amplified from vacuum fluctuations, as its wavelength 
increases from sub-Hubble to super-Hubble scales. After in- 
flation ends, the universe reheats and leads to the production 
of charged particles leading to a dramatic increase in the 
plasma conductivity. Then the electric field E would get 
shorted out while the magnetic field B of the EM wave gets 
frozen in. This is the qualitative picture for the generation 
of primordial fields during the inflationary era. 

There is however another potential difficulty; since a(t) 
is almost exponentially increasing during slow roll infla- 
tion, the predicted field amplitude, which behaves say as 
B ~ 1 /a e is exponentially sensitive to any changes of the 
parameters of the model which affects e. Therefore models 
of magnetic field generation can lead to fields as large as 
B ~ 10~ 9 G (as redshifted to the present epoch) down to 
fields which are much smaller than that required for even 
seeding the galactic dynamo. For example in model consid- 
ered by Ratra (1992) with f 2 (4>) ~ e"^, with (j) being the 
inflaton, one gets B ~ KT 9 to 1(T 65 G, for a ~ 20 - 0. 
Note that the amplitude of scalar perturbations generated 
during inflation is also dependent on the parameters of the 
theory and has to be fixed by hand. But the sensitivity to 
parameters seems to be stronger for magnetic fields than for 
scalar perturbations due to the above reason. Nevertheless 



2 Interestingly slower decay can occur for modes comparable to the cur- 
vature scale, in open FRW models due to coupling to the spatial curvature 
(Barrow and Tsagas, 2008). 



one may hope that there would arise a theory where the pa- 
rameters are naturally such as to produce interesting primor- 
dial magnetic field strengths. We describe below one frame- 
work for magnetic field generation during inflation, keeping 
the discussion quite general without specifying any specific 
inflation model. A nice treatment of inflationary generation 
of magnetic fields, which we follow to some extent, is given 
by Martin and Yokoyama (2008), where some specific mod- 
els are also discussed. 

3.1 Quantizing the EM field 

Let us assume that the scalar field (j> in Eq. (1431 is the field 
responsible for inflation and also assume that this is the sole 
term which breaks the conformal invariance of the electro- 
magnetic action. The total action is given by 

1 



s = -— i d 4 Iv q[ s a Y7 2 (^)^^] 
i 



16tt 

d 4 x 



-g^d^d^ + V(<t>) 

Maxwell equation now become [f 2 F^ u ]^ = 0, or 

1 d 

^gdx v 

The scalar field satisfies 



(44) 



(45) 



1 



d 



dV 



]_d[_ 
8tt dxj> 



V-9 dx" 

We assume that the electromagnetic field is a 'test' field 
which does not perturb either the scalar field evolution or 
the evolution of the background FRW universe. We take the 
metric to be spatially flat with 

ds 2 = -dt 2 + a 2 [dx 2 + dy 2 + dz 2 ] 

= a 2 {rj) [-dr? + dx 2 + dy 2 + dz 2 ] (47) 

where r\ = j (dt/a) is the conformal time. Throughout the 
discussion in this section, we use conformal time (rf) and 
co-moving space (x, y, z) as our four co-ordinates and all 
tensors when explicitly specified will be in this co-ordinate 
frame. It is convenient to adopt the Coulomb gauge: 

Ao(7?,x)=0, d j A j (r l ,x)=0. (48) 

In this case the time component of Eq. d43T > becomes a trivial 
identity, while the space components give 

J' 



A" 



2-4 



a 2 d^Ai 







(49) 



where we have defined — dk = a -2 ^ 7 <9ft, and a 
prime denotes derivative with respect to 77. In fact a 2 djd^ is 
the usual spatial V 2 operator with respect to the co-moving 
spatial co-ordinates. 

We can also use Eq. © to write the electric and mag- 
netic fields in terms of the vector potential. Note that the 
four velocity of the fundamental observers used to define 
these fields is now given by w M = (1/a, 0, 0, 0). The time 
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components of and B^ are zero, while the spatial com- 
ponents are given by 

Ei = --4; Bi = -e* jk P l 6 km (diA m ) (50) 
a a J 

Note that these spatial components are the same as when 
t is used instead of 77 as the time co-ordinate. This is be- 
cause the transformation t — >• 77 is independent of x, and we 
have not transformed the space co-ordinates. For a constant 
/ Eq. ( |49l shows that A, simply satisfies the usual wave 
equation in r\ and x co-ordinates, whose solutions are plane 
waves with constant amplitude. Then the amplitude of Bi 
then scales as 1 /a, while that of B 1 scales as 1 /a 3 and so 
the amplitude of B a scales as 1 /a 2 as before. 

We would like quantize the electromagnetic field in the 
FRW background. For this we treat Aj as the co-ordinate, 
and find the conjugate momentum II 1 , by varying the EM 
Lagrangian density Cem = —pF^F^KlQn), with re- 
spect to A! i . We get 

5Cem 



rr 



SAL 



in 



fa 2 g^A', IT 



in 



A a 4 



To quantize the electromagnetic field, we promote A 1 and 
Ilj to operators and impose the canonical commutation re- 
lation between them, 

= iS i Xj (x-y). (51) 

Here the term Pj(fc) = {5) - S om {k l k m /k 2 )) is introduced 
to ensure that the Coulomb gauge condition is satisfied and 
S' l ± . is the transverse delta function. This quantization con- 
dition is most simply incorporated in Fourier space. We ex- 
pand the vector potential in terms of creation and annihila- 
tion operators, b\(k) and b\(k), with k the co-moving wave 
vector, 



A 1 (03,77) = V 47T 



d 3 k 



(2tt 



E e A( fc ) 



b x {k)A(k,T])e 



k x 



b[(k)A*(k, v )e~ 



ikx 



■(52) 



Here the index A = 1, 2 and e\(k) are the polarization vec- 
tors, which form part of an orthonormal set of basis four- 
vectors, 



a I \ a 



,e£= 0, 



k 



(53) 



The 3-vectors e\ are unit vectors, orthogonal to k and to 
each other. The expansion in terms of the polarization vec- 
tors incorporates the Coulomb gauge condition in Fourier 
space. It also shows that the free electromagnetic field has 
two polarization degrees of freedom. If we substitute the 
Fourier expansion given in Eq. ( l52l into Eq. d49b . we find 
that the Fourier coefficients A — (aA(k, 77)) satisfy, 



f - 

A" + 2jA' + k 2 A = 



(54) 



One can also define a new variable A = a(rj)f(jj)A(r],k) 
in order to eliminate the first derivative term, to get 



A"( V ,k)+[ k z - 



r 
f 



A(v,k) = 



(55) 



We can see that the mode function A satisfies the equa- 
tion of a harmonic oscillator with a time dependent mass 
term. The case / = 1 corresponds to the standard EM ac- 
tion where A oscillates harmonically with time. One needs 
/" > to have possible growth of magnetic fields. 

The quantization condition given in Eq. dBTT i is imple- 
mented by imposing the following commutation relations 
between the creation and annihilation operators, 



bx(k),b{,(k')] = (2n) 3 5 3 (k-k')S X y 
[b x (k),b x ,(k')] = \bl(k),bl,(k')] = 0. 



(56) 



We also define the vacuum state |0 > as one which is anni- 
hilated by b x (k), that is b\(k)\0 >= 0. Note that the choice 
of the initial quantum state will be decided by the choice of 
the mode function A as below. To check if Eq. ( BTT l is indeed 
satisfied, we can substitute the Fourier expansion of A 1 and 
into the commutator [A 1 , II,-], and use the commutation 
relations Eq. d56l ). We get 

7^3 £4(*)ejA(fc) x 

^ ' A=l 

e ik < x -y^W(k, V )fa 2 , (57) 
where we have defined the complex Wronskian, 

W{k,r]) = [A A'* - A* A']. 

One can check that the polarization four-vectors satisfy the 
completeness relation, 



2 

J2e\(k)e jX (k) = P;(k) 
\=i 



(58) 



Let us also define the Wronskian associated with A, given 
by W = [A A'* - A* A'}. We have W = a 2 W, and since 
A satisfies Eq. (O, we g et W = -(2f / f)W. Integrating 
this equation we get W — a 2 W cx (l// 2 ). Substituting 
the expression for W into Eq. d57l >. and using Eq. ( fSSb . we 
can verify that the quantization condition Eq. (IBTI ) is indeed 
satisfied, provided we fix the constant of proportionality in 
W such that W = {i/f 2 a 2 ). 

Once we have set up the quantization of the EM field, it 
is of interest to ask how the energy density of the EM field 
evolves. The energy momentum tensor is given by varying 
the EM Lagrangian density with respect to the metric, 



T, 



2 S[^C 



EMI 



flu 



in 



(59) 
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The energy density T^u^u" can be written as the sum of a 
magnetic and electric contributions. The energy density due 
to the magnetic part of the EM field is given by 

P 



16tt 



[A m ,% - A^ m ][Aij - Ajj]g ij g ml 



i B;B l 



(60) 



Similarly, the energy density due to the electric field is 

T>»u» = ^A'W = . (61) 

We substitute the Fourier expansion of Ai into Eq. (l60b and 
Eq. drjTT ). and take the expectation value in the vacuum state 
|0 >. Let us define 

p B =<0\T^u^\0>, p E =< 0\T^u u \0 > . 

Using the properties 

b x (k)\0> = 0, 

<0|6 A (fe)6|,(p)|0> = (2ir) 3 6(k - p)6xx>, (62) 

we get for the spectral energy densities in the magnetic and 
electric fields, 



dpB 
<iln k 

dps 
dink 



1 

2^ 



f 2 k 3 
2tt 2 a 4 



k\A(k, V )f 



A(k,r]) 
f 



(63) 



(64) 



Once we have calculated the evolution of the mode function 
A(k,rf), the evolution of energy densities in the magnetic 
and electric parts of the EM field can be calculated. 

3.2 Evolution of normal modes 

Consider for example a case where the scale factor a{rj) 
evolves with conformal time as 

, 1+0 



a(rj) = a 



(65) 



The case when f3 = — 2 corresponds to de Sitter space- 
time of exponential expansion in cosmic time, or a(t) oc 
exp (Ht). On the other hand for an accelerated power law 
expansion with a(t) = a^(t/t^) p and p > 1, integrating 
dt = adrj, we have 



to 



-i/(p-i) 



a{ri) 



«o 



-p/Cp-i) 



(66) 



a (p- 1) 

'-o (p - l)r/" 
to 

Here we have assumed that r\ — > 0_ as t — > oo, such that 
during inflation, the conformal time lies in the range — oo < 
r\ < 0. In the limit of p ^> 1, one goes over to an almost 
exponential expansion with (3 — >• — 2 — 1/p. 

Let us also consider a model potential where the gauge 
coupling function / evolves as a power law, 

f( V ) ex a a . (67) 



This could obtain for example for exponential form of f(<fi) 
and power law inflation. We then have 

f" _ 7(7-1) 
/ 



7 = 0(1 + 0) 



(68) 



Then the evolution of the mode function A is given by 

.2 7(7 - 1) 



A"(k,ri) 



Tj 



A{k, V )=0, 



(69) 



whose solution can be written interms of Bessel functions, 

A = (-k v ) 1 / 2 [ d(fc)V 1/2 (-fc,) 

+ C 2 (k)J_ 1+1/2 (-kr))], (70) 

where Ci(fc) and C2(fc) are scale-dependent coefficients to 
be fixed by the initial conditions. 

Let us define the Hubble radius at any epoch as the 
length scale Rh — 1 /H (in units where the speed of light 
c = 1). The initial conditions to determine the constants in 
Eq. ( T70b are specified for each mode (or wavenumber fc), 
when it is deep within the Hubble radius. That is when the 
proper scale length associated with the mode (k/a)^ 1 is 
much smaller than Rh, or (k/aH) 3> 1. For such small 
scale modes one assumes that the effects of space-time cur- 
vature are negligible and thus the mode function goes over 
to that relevant for the Minkowski space vacuum. Recall that 
the expansion rate is given by H(t) = a /a = a' /a 2 . (Here 
and below A is derivative of A with respect to proper time.) 
For the expansion factor given in Eq. d66t . we have a' /a = 
-(p/(p - 1))(1/t/), and forp > 1, aH -> -I/77. Thus the 
ratio the Hubble radius to the proper scale of a perturbation 
is given by (l/i?)(a/fc) _1 = k/(aH) = —kr\. A given 
mode is therefore within the Hubble radius for —kr\>\ 
and outside the Hubble radius when — krj < 1. 

In the short wavelength limit, (k/a)/H — (—krj) — > 
00, the solutions of Eq. d69l are simply A oc exp (±ikrj). 
The assumption that the gauge field for these modes is in the 
Minkowski space vacuum state leads us to pick the 'pos- 
itive' frequency mode A = cq exp (— ikrf), and from the 
Wronskian condition that W = i/f 2 a 2 , the constant cq is 
fixed to Co = l/y2k. Thus we assume as initial condition 
that as (—krj) — > 00, 



A 



1 



exp 



— ikrj 



(71) 



This fixes the constants in Eq. d70t to be, 

7r exp (—inj/2) 



cm 



C 2 (k) 



Ak cos(7T7) 

7r exp (z7r(7 + l)/2) 



4fc cos(7T7) 

where we have used the asymptotic expansion, that for x 

00, 



(72) 



J v {x) -> 



(v + -)- 



In the opposite limit of modes well outside the Hubble ra- 
dius, or at late times, with {—k-q) — > 0, we get from Eq. d70l >, 

A -+ k- 1 / 2 [ Cl ( 7 )(-M 7 + c 2 ( 7 )(-fcr;) 1_7 ] , (73) 
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where 



ci 



27+1/2 r( 7 + |)cos(7r 7 ) : 



C-2 



23/2-7 T(| - 7 ) C0S(^ 7 ) : 

Here we have used the property that 



(74) 



J v (x) 



as x — >• 



2 v r(i/ + i)- 

From Eq. ( f73b one sees that the ci term dominates for 7 < 
1/2, while C2 term dominates for 7 > 1/2. 

As an aside, we note that the late time (small k) solution 
can also be got for a more general / from directly integrat- 
ing Eq. d55l > in the limit k — > 0. We get 

drj 

T 2 

If we substitute / oc a a and a cx into Eq. ( l75l ). we 

will recover the long wavelength solution given in Eq. ( 1731 ). 



A->c\f + c 2 f , _> . 



(75) 



3.3 The generated magnetic and electric fields 

We can now calculate the spectrum of ps and pe in the late 
time, super Hubble limit. Substituting Eq. ( f73b into Eq. d63K 
we get for the magnetic spectrum, 

dp B _ T(n) 4 ' ! - x 



din A; 



H' 

2tt 2 V aH 



2tt 2 



H 4 (-A77) 



4+2n 



(76) 



where rt = 7if7<l/2 and n = 1 — 7 for 7 > 1/2, and 

7T 



J-(n) 



(77) 



2 2 »+ 1 r 2 (n+|)cos 2 (7rn) 

During slow roll inflation, the Hubble parameter H is ex- 
pected to vary very slowly, and thus most of the evolution of 
the magnetic spectrum is due to the (—krf) A+2n factor. One 
can see that the property of scale invariance of the spectrum 
(with 4 + 2n = 0), and having ps ~ a go together, and 
they require either 7 — 3 or 7 = — 2. 

We can also calculate the electric field spectrum in a 
very similar manner. We first find {A/f)' from Eq. d70b . 
using the identities 



V _ 

*' v ,J V 
X 



X 



and then take the limit (—kf]) —> 0. Substituting the result 
into Eq. d64l gives 

4+2m 



dps _ Q{m) H ± 
dink 2tt 2 \aH 



£( TO ) i/4 1 1 ^4+2m 



2tt 2 



(78) 



where now m = 7 + 1 if 7 < —1/2 and m = —7 for 
7 > -1/2, and 

7T 



Q(m) 



2 2m + 3 r 2 (m + f ) cos 2 (7rm) 



(79) 



Thus having a scale invariant magnetic spectrum implies 
that the electric spectrum is not scale invariant, and in addi- 
tion can vary strongly with time. For example if 7 = 3, then 
(4 + 2m) = -2, although (4 + 2n) = 0. In this case, at late 
times as {—k-q) — > 0, the electric field increases rapidly, 
with pe oc {—kri)~ 2 — > 00. There is then the danger of 
its energy density exceeding the energy density in the uni- 
verse during inflation itself, unless the scale of inflation (or 
the value of H A is sufficiently small. Such values of 7 are 
strongly constrained by the back reaction on the background 
expansion they imply (Martin and Yokoyama, 2008). 

On the other hand consider the case near 7 = —2. In 
this case the magnetic spectrum is scale invariant, and at the 
same time (4 + 2m) = 2, and so the electric energy density 
goes as {—krj) 2 —> as (— krj) — > 0. Thus these values of 7 
are acceptable for magnetic field generation without severe 
back reaction effects. 

We now discuss the evolution of the field after inflation. 
Post inflationary reheating is expected to convert the energy 
in the inflaton field to radiation (which will include vari- 
ous species of relativistic charged particles). For simplic- 
ity let us assume this reheating to be instantaneous. After 
the universe becomes radiation dominated its conductivity 
(a) becomes important. Indeed Turner and Widrow (1988) 
showed that the ratio a/H 3> 1. In order to take into ac- 
count this conductivity, one has to reinstate the interaction 
term in the EM action, given in Eq. (JTJ. Further, as the in- 
flaton has decayed, we can take / to have become constant 
with time and settled to some value Jq. Varying the action 
with respect to now gives 

4nJ^ 



PC" 

• 11 



/o 2 



The value of fa thus goes to renormalize the value of electric 
charge e to be e^v = e/f§- This aspect raises an additional 
potential problem for values of 7 w —2, which has been re- 
cently emphasized by Demozzi, Mukhanov and Rubinstein 
(2009) (DMR). 

Suppose the inflationary expansion is almost exponen- 
tial with f3 = —2, then for 7 —2, we have a = 7/(1 + 
j3) 2. This implies that the function / = f i (a/ai) 2 in- 
creases greatly during inflation, from its initial value of fi 
at a = <Zj. Thus if we want /q ~ 1 at the end of inflation, 
then at early times /, <C fo and the renormalized charge at 
these early times e^v = e/ f 2 3> e. DMR argue that one is 
then in a strongly coupled regime at the beginning of infla- 
tion where such a theory is not trustable. There is however 
the following naive caveat to the above argument: Suppose 
one started with a weakly coupled theory where fi ~ 1. 
Then at the end of inflation fo ^> /j, and so the renormal- 
ized charge e^r <C e. Such a situation does not seem to 
have the problem of strong coupling raised by DMR; how- 
ever it does leave the gauge field extremely weakly coupled 
to the charges at the end of inflation. This also means that 
even if ps is large, the magnetic field strength itself as de- 
duced from Eq. d60l is B l Bi = 8ttpb/ fo <C 8npB- Ideas 
to sort out this difficulty need to be explored, whether for 
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example one can relax a large fo back to /j, without now 
re-generating strong electric fields. 

Let us proceed by assuming that we have absorbed /q 
into e. In the conducting plasma which obtains after reheat- 
ing, the current density will be given by the Ohm's law of 
Eq. d27b . The fluid velocity at this stage is expected to be 



that of the fundamental observers, i.e. 



Thus the 



spatial components J 1 = aE 1 = —g^Aj. Let us assume 
that the net charge density is negligible and thus neglect 
gradients in the scalar potential A$. Then the evolution of 
the spatial components of the vector potential is given by 

Ai + (H + 47rcr)i j; - dj&Ai = . (80) 

We see that any time dependence in Ai is damped out on 
the inverse conductivity time-scale. To see this explicitly, 
consider modes which have been amplified during inflation 
and hence have super Hubble scales k/(aH) <C 1. Also let 
us look at the high conductivity limit of cr/H ^> 1. Then 
Eq. d80l reduces to 

Ai + 4%aAi = 



whose solution is given by 

Di(as 



Ai 



-4TT<7t 



47TCT 



+ D 2 (x) 



(81) 



We see that the D\ term decays exponentially on a time- 
scale of (47rcr) _1 -c (1/H). This leaves behind a constant 
(in time) Ai = D2(x). Thus the electric field Ei = 0, and 
the high conductivity of the plasma has led to the shorting 
out of the electric field. Note that the time scale in which the 
electric field decays does not depend on the scale of the per- 
turbation, that is the a dependent damping term in Eq. d80l 
has no dependence on spatial derivatives. As far as the mag- 
netic field is concerned, Eq. d50l shows that Bi ~ 1/a when 
Ai = D2(x). Therefore Bi ~ 1/a 2 , as expected when the 
magnetic field is frozen into the highly conducting plasma. 

Let us now make a numerical estimate of the strength 
of the magnetic fields generated in the scale invariant case. 
For both 7 = — 2 and 7 = 3, we have from Eq. d76l l and 
Eq. {77), 

dp B 9 



(82) 



din A; An 2 ' 

Cosmic Microwave Background limits on the amplitude of 
scalar perturbations generated during inflation, give an up- 
per limit on H/M v i ~ 10~ 5 (cf. Bassett, Tsugikawa and 
Wands, 2006). Here M pl = 1/y/G is the Planck mass. The 
magnetic energy density decreases with expansion as 1/ a 4 , 
and so its present day value ps(0) = Ps(a/ /ao) 4 , where 
ay is the scale factor at end of inflation, while ao is its 
present day value. Let us assume that the universe transited 
to radiation domination immediately after inflation and use 
entropy conservation, that is the constancy of gT 3 a 3 during 
its evolution, where g is the effective relativistic degrees of 
freedom and T the temperature of the relativistic fluid. We 
get 

1/4 



ao 
a f 



gift! ff l/ 2M V 2 



1/3 

9o 



T 



90 

8^ 



Taking g f ~ 100, gives (a /a f ) ~ 10 29 (iJ/10- 5 Af p; ) 1/2 > 
This leads to an estimate the present day value of the mag- 
netic field strength, Bq at any scale, 

B Q ~ 5 x 10- 10 G I J ) . (83) 

Thus interesting field strengths can in principle be created if 
the parameters of the coupling function / are set appropri- 
ately and the problems highlighted by DMR can be circum- 
vented. Note that the strength of the generated field is sensi- 
tive to even slight departures from scale invariance. Suppose 
7 = -2 + e, with e < 1, then 

2e 



9 



dp B 



d In k 47r 2 



aH 



(84) 



valid for (k/aH) <C 1. Assuming a radiation dominated 
universe immediately after inflation, and matter domination 



from a redshift z P 



(k/aH) - 3 x 10 



3300, we estimate 
_24 ' ^ 



H 



1/2 



,1/iMpcTV \lO~ 5 M pli 

Thus at galactic scales of k = lh Mpc -1 , B will be smaller 
or larger by a factor of ~ 10~ 5 , if one takes e = ±0.2 re- 
spectively. This shows the sensitivity of the magnetic field 
amplitude to small changes in the parameters of any gener- 
ation model, as mentioned earlier. 

4 Discussion 

We end with a few comments. We have emphasized the use 
of tetrads in defining properly behaved magnetic and elec- 
tric fields. Such a procedure is also followed when using 
with the 3 + 1 formalism in the context of black-hole elec- 
trodynamics (cf. Macdonald et al., 1986). Thus it seems 
to us somewhat surprising why this does not usually get 
mentioned when dealing with MHD in the context of cos- 
mology. Regarding magnetic field generation during infla- 
tion, we have a paradigm, but no compelling model as yet. 
Clearly more work is required to find such a model, which 
at the same time avoids the problems of back reaction, the 
strong coupling problem mentioned by DMR, or the strong 
decay of coupling constants. 

Another possibility is generation of primordial fields is 
during a later phase transition, like the Electro-Weak or the 
quark-hadron transition. Here the main problem is that the 
generated field typically has a tiny correlation scale, less 
than the Hubble radius H~ l , at the epoch of the phase tran- 
sition, unless magnetic helicity is also generated. Interest- 
ingly there are several ideas for such helicity generation dur- 
ing the Electro- weak phase transition (cf. Vachaspati, 2001; 
Diaz-Gil et al., 2008; Copi et al., 2008). Magnetic energy 
decay conserving helcity can then lead to an inverse cas- 
cade and larger coherence scales, as first emphasized in the 
early universe context by Brandenburg, Enquist and Ole- 
sen (1996) (see also Christensson, Hindmarsh and Branden- 
burg, 2001; Banerjee and Jedamjik, 2005), but that is a story 
for another review. 
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